home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 1998 November / IRIX 6.5.2 Base Documentation November 1998.img / usr / share / catman / p_man / cat3 / complib / CPPDI.z / CPPDI
Text File  |  1998-10-30  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. CCCCPPPPPPPPDDDDIIII((((3333FFFF))))                                                            CCCCPPPPPPPPDDDDIIII((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CPPDI   - CPPDI computes the determinant and inverse of a complex
  10.      Hermitian positive definite matrix using the factors computed by CPPCO or
  11.      CPPFA .
  12.  
  13.  
  14. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  15.       SUBROUTINE CPPDI(AP,N,DET,JOB)
  16.  
  17. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  18.      On Entry
  19.  
  20.      AAAAPPPP COMPLEX (N*(N+1)/2)
  21.         the output from CPPCO or CPPFA.
  22.  
  23.      NNNN INTEGER
  24.         the order of the matrix  A .
  25.  
  26.      JJJJOOOOBBBB INTEGER
  27.         = 11   both determinant and inverse.
  28.         = 01   inverse only.
  29.         = 10   determinant only.  On Return
  30.  
  31.      AAAAPPPP the upper triangular half of the inverse .
  32.         The strict lower triangle is unaltered.
  33.  
  34.      DDDDEEEETTTT REAL(2)
  35.         determinant of original matrix if requested.
  36.         Otherwise not referenced.
  37.         Determinant = DET(1) * 10.0**DET(2)
  38.         with  1.0 .LE. DET(1) .LT. 10.0
  39.         or  DET(1) .EQ. 0.0 .  Error Condition
  40.  
  41.      AAAA division by zero will occur if the input factor contains a zero on the
  42.      diagonal and the inverse is requested.  It will not occur if the
  43.      subroutines are called correctly and if CPOCO or CPOFA has set INFO .EQ.
  44.      0 .  LINPACK.  This version dated 08/14/78 .  Cleve Moler, University of
  45.      New Mexico, Argonne National Lab.  Subroutines and Functions BLAS
  46.      CAXPY,CSCAL Fortran CONJG,MOD,REAL
  47.  
  48.  
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.